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Abstract 

A detailed next-to-next-to- leading order (NNLO) QCD analysis is performed for 
the experimental data of the CCFR collaboration for the xF^ structure function. The- 
oretical ambiguities of the results of our NNLO fits are estimated by application of 
the Fade resummation technique and variation of the factorization and renormal- 
ization scales. The NNLO and N^LO as{Q^) M5-matching conditions are used. 
In the process of the fits we are taking account of twist-4 1/Q^-terms. We found 
that the amplitude of the a;-shape of the twist-4 factor is decreasing in NLO and 
NNLO, though some remaining twist-4 structure seems to retain in NNLO in the 
case when only statistical uncertainties are taken into account. The question of the 
stability of these results to the application of the [0/2] Fade resummation technique 
is considered. Our NNLO results for as{Mz) values, extracted from the CCFR xF^ 
data, are a^iMz) = 0.118 ± 0.002(siai) ± 0m5{syst) ± 0.003{theory) provided the 
twist-4 contributions are fixed through the infrared renormalon model and as{Mz) — 
O.Ulto^oilistat) ± 0.005(syst) ± 0.003{theory) provided the twist-4 terms are consid- 
ered as free parameters. 
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1. Introduction. 

Deep-inelastic lepton-nucleon scattering (DIS) belongs to classical and constantly stud- 
ied processes in modern particle physics. The traditionally measurable characteristics of 
vN DIS are SFs (structure functions) F2 and xF^. It should be stressed that the program 
of getting information about the behavior of SFs of vN DIS is among the aims of the 
experimental program of Fermilab Tevatron and CCFR/NuTeV collaboration, in partic- 
ular. The CCFRR collaboration started to study the vN scattering process in 1980 Q. 
The data for SFs of fN DIS, obtained by a follower of the CCFRR collaboration, namely, 
the CCFR group, were distributed among the potential users in the beginning of 1997 ||2|, 
while the final results of the original CCFR DGLAP NLO analysis of these data were 
presented in the journal publication of Ref.|H]. 

This experimental information was already used in the process of different NLO anal- 
yses, performed by CTEQ, MRST, and GRV groups (see Refs.Q, ^, ^, respectively). 
Subsequent steps of performing NLO and the first NNLO analysis of the CCFR data with 
the help of the Jacobi polynomial - Mellin moments version of the DGLAP method were 
made in Refs.||7|- ||l^ (definite stages in the development of this formalism are described 
in Refs.ll-ll]). 

In the process of the analysis of Refs.[P]-[12|, the authors used important information 
about the NNLO corrections to the coefficient functions |16] and results of the complicated 
analytic calculations of the NNLO corrections to the anomalous dimensions of the NS 
moments with n = 2, 4, 6, 8, 10 [17|, supplemented with the estimated values of the NNLO 
coefficients of the anomalous dimensions of n = 3, 5, 7, 9 NS moments, which were obtained 
in Ref.[^ with the help of the smooth interpolation procedure, previously proposed in 
Ref . [Il8| . Moreover, the attempts to obtain the shape of the twist-4 contributions and to 



study the predictive abilities of the IR-renormalon (IRR) model of Ref.|l£] were made 
(for certain details of modeling the effects of the power-suppressed contributions to XF3 
and other measurable physical quantities in the IRR language see Ref. 1 20] and Ref.|21| 
respectively, for more details, see the review of Ref. ]22|j ) . 

However, an important question of estimating theoretical uncertainties of the NNLO 
analysis of the CCFR data of Ref.||ll| was not analyzed in detail. These uncertainties can 
be specified in part after taking into account 

1) differences in the definitions of as{Q'^) matching conditions (see e.g. Refs.[23, 24, 
pSj [2^), which are responsible for penetrating into the energy region, characteristic of 
/ = 5 numbers of flavours, where the pole of the Z'^-boson manifests itself; 

2) NNLO QCD contribution to the matching condition of Ref. ||23[| (corrected in Ref. [^] ) , 



the 4-loop coefficient of the QCD /3-function |28|, which is entering into the N^LO approxi- 
mation of the renormalization group evolution equation for the Mellin moments, the N^LO 
expression for the coupling constant and to the calculated N'^LO-term ||29|| of the matching 
condition; 

3) theoretical uncertainties due to non-calculated N^LO contributions to the coefficient 
functions and the anomalous-dimension functions; 

4) in the NNLO it is also desirable to carefully analyse the dependence of the results 
obtained on the choice of the Jacobi polynomial parameters a and (3 (see definitions 
below). These parameters are entering into the theoretical expression for the reconstructed 



2 



structure function (in the NLO the problem of their fixation was studied in detail in 
Ref ■ []l5| ) . This analysis is of relevant importance in view of doubts in the applicability of 
the Jacobi evolution method (see Ref. pC|] ) , which, however, were immediately dispelled in 
Ref.P; 

5) last, but not least, uncertainties are related to problems of the sensitivity of the 
obtained results to the choice of factorizatation and renormalization scales. It is worth 
noting that these problems are in relation with fixing the ambiguities due to already 
calculated and uncalculated N^LO QCD effects, which can be modelled using the Fade 
resummation technique. 

This work is devoted to the analysis of important problems outlined above and to a 
more detailed determination of values of as{Mz) and the x-shape of the twist-4 power- 
suppressed term in available orders of perturbative QCD taking into account the effects 
listed above. We supplement the NNLO fits of Ref. [|lT| by the N'^LO analysis based on the 
application of the Fade resummation technique (for a review, see Ref. [^] ) , developed in 
QCD in a definite form in Refs.jS^, U], and considered previously as a possible method of 
fixing theoretical uncertainties in the analysis of DIS data in Ref.||3^. It should be stressed 
that a posteriori this technique gives results similar to those obtained by different methods 
of fixing scale-scheme dependence ambiguities (compare the results of Ref. [^] with the 



results of Refs.|33, ^ obtained by using the Fade resummation technique). Thus, our 
analysis could be considered as an attempt to estimate perturbative QCD uncertainties 
beyond the NNLO level. Moreover, it could give us a hint whether the outcomes of 
our NNLO fits, related to perturbative and non-perturbative sectors, are affected upon 
including the explicitly calculated and estimated N^LO QCD effects. 

Another new important ingredient of our analysis, discussed in brief in Ref.[^], is the 
analysis of the problem of the initial scale choice. In particular, we will demonstrate that 
due to unnaturally large NNLO corrections to the renormalization-group improved n=2 
Mellin moment of xF^, , it is essential to choose the value of the initial point in a vicinity 
of the scale Ql=2Q GeV^ . 

2. The theoretical background of the QCD analysis. 

Let us define the Mellin moments for NS SF xF^{x, Q^): 

= /' x"-^F3(x, Q^)dx (1) 
Jo 

where n = 2,3,4, .... These moments obey the following renormalization group equation 

^^§-^ + Pi^s)^^+l'^;:l{As))M^^{Q'/^^\As{^?)) = Q (2) 

where Ag = as/(47r). The renormalization group functions are defined as 

/^^^ = ^nI{A,) = E l%^n)Ar (3) 



j>0 
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where Z^^ are renormalization constants of the corresponding NS operators. The solution 
of the renormahzation group equation can be represented in the following form : 



exp 



dx 



(4) 



where M^^(Qq) is a phenomenological quantity related to the initial- scale dependent 
factor. At a fixed momentum transfer Qq, it will be parameterized in the simple form 



(5) 



with 7 7^ or 7 = 0. It is identical to the form used by the CCFR collaboration Q. 
In principle, following the models of parton distributions used in Refs.[^, |6[, one can in 
Eq.(5) add a term proportional to y/x. However, since this term is important only in the 
region of rather small x, we neglect it in our analysis. 

In N^LO, the coefficient function can be defined as 

cPsiAs) = 1 + C7(i)(n)^ + C(2)(n)^2 ^ c(3)(n)A3, 

while the corresponding expansion of the anomalous-dimension term is 



exp 



P{x) 



dx 



where 

AD{n,As) = [l+p{n)A,{Q^)+q{n)A,{Q''f + r{n)As{Q^f] 
and p{n), q{n) and r{n) have the following form: 



2\ Pi 

k 

(3o 



,(n) = i ( 2pinr + + ^P,in)'J^i^^ - j^Un)^ 



/3o JPo 
(/3?-/32/3o) 



r(n) 



2p(n)^ -|- 6p{n)q(n) + 



/3o 



Pi 



(6) 

(7) 
(8) 

(9) 

(10) 

(11) 



Pi 



Po Po Po Po Ho 

The coupling constant As{Q'^) can be expressed in terms of the inverse powers of L = 
ln(QVA|^) as A^^o ^ j^LO + A^f ^ j^^nlo + ^a^^nlo j^fLO ^ 

A^^^^o ^ ^j^fLO ^ where 

^ (12) 



A^o 



PoL 
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NNLO 



(13) 



N^LO 



-±-^[Pf(-lnHL) + ^ln\L) + 2/n(L) - 1) 
-3/3o/3i/32/n(L)+/32^] . 



(14) 



Notice that in our normalization the numerical expressions for /3o, Pi, P2 and P3 read 



/So 
Pi 
P2 
Pz 



11 - 0.6667/ 
102 - 12.6667/ 

1428.50 - 279.611/ + 6.01852/^ 

29243.0 - 6946.30/ + 405.089/^ + 1.49931/' 



(15) 



where the expression for P2, was obtained in Ref.|28|. The inverse- log expansion for 
Avl^^^*-^ which incorporates the information about the coefficient P^ was presented in 
Ref.O. 

A few words are to be said about the approximation used for the anomalous-dimension 
function 7^g(As). The analytic expression for its one- loop coefficient is well-known: 



(8/3)[4S^=i(l/j) - 2/n(n + 1) - 3]. In the cases of both F2 and xFg, nu- 



merical expressions for 7]y^(n)-coefficients are given in Table 1. 



n 


77VS,F2(^) 


7!v5,F3(^) 




7ivs(^)l[i/i] 


7jvs(^)l[o/2] 


2 


71.374 


71.241 


612.1 


5259 


5114 


3 


100.801 


100.782 


837.4 


6959 


6900 


4 


120.145 


120.140 


1005.8 


8421 


8414 


5 


134.905 


134.903 


1135.8 


9563 


9562 


6 


147.003 


147.002 


1242.0 


10493 


10482 


7 


157.332 


157.332 


1334.0 


11310 


11280 


8 


166.386 


166.386 


1417.5 


12077 


12012 


9 


174.468 


174.468 


1493.5 


12784 


12706 


10 


181.781 


181.781 


1559.0 


13370 


13271 


11 


188.466 


188.466 


? 


? 


? 


12 


194.629 


194.629 


? 


? 


? 


13 


200.350 


200.350 


? 


? 


? 


14 


205.689 


205.689 


? 


? 


? 



Table 1. The used numerical expressions for NLO and NNLO coefficients of anomalous 
dimensions of moments of the NS SFs at / = 4 number of flavours and the N'^LO Fade 
estimates. 



These results are normalized to the case of / = 4 numbers of active flavours. In the 
same Table we present the numerical expressions for 7]vs(^)) used in the process of fits. 
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(2) 

In the cases of n = 2,4,6,8,10 they follow from the explicit calculations of 7]v5F2(^)" 
terms [p!7|, normalized to / = 4, while the n = 3,5,7,9 numbers are fixed by using the 



smooth interpolation procedure originally proposed in Ref.|18]. Note in advance, that since 

(2) (2) — 

INS Fg (f^)-coefficients differ from jj^g p (n)-terms, though by presumably small additional 
contributions (for discussions, see Ref.pl), it would be interesting to verify the accuracy of 
the expression for used in the process of our NNLO 2:^3 fits, by explicit analytic 

calculations of the NNLO contributions to anomalous dimensions of odd moments of the 
structure function. 

Let us now describe the procedure of fixing other theoretical uncertainties. After the 
work of Ref.|33] it becomes rather popular to model the effects of higher order terms of 
perturbative series in QCD by the expanded Fade approximants. 

In the framework of this technique, the values of terms C^^\n) and r(n) could be 
expressed as 



Fade [1/1] : C^'^Xn) = [C^'^\n)f / C^^\n) (16) 

r(n) = qin)"^ /p{n) (17) 

Fade [0/2] : C'-^\n) = 2C''^\n)C^^\n) - [C'-^Hn)]^ (18) 

r(n) = 2p{n)q{n) - [pin)f (19) 

The numerical values for p{n) and q{n), obtained from the results of Table 1 and 
definitions of Eqs.(9)-(11), together with the values of the coefficients C^^\n) and C^'^\n) 
(which come from the calculations of Ref . []l6| ) , are presented in Table 2. 



n 


p{n) 


q{n) 


r{n)\[i/i] 


r{n)\[o/2] 


CW(n) 


C7i^)(n) 




C^'Kn)\[o/2] 


2 


1.646 


4.232 


10.829 


9.476 


-1.778 


-47.472 


-1268 


174 


3 


1.941 


4.774 


11.738 


11.218 


1.667 


-12.715 


97 


-47 


4 


2.050 


5.546 


15.003 


14.123 


4.867 


37.117 


283 


246 


5 


2.115 


6.134 


17.790 


16.486 


7.748 


95.408 


1175 


1013 


6 


2.165 


6.595 


20.087 


18.407 


10.351 


158.291 


2421 


2168 


7 


2.210 


7.039 


22.421 


20.318 


12.722 


223.898 


3940 


3638 


8 


2.252 


7.525 


25.138 


22.471 


14.900 


290.884 


5679 


5360 


9 


2.294 


8.018 


28.027 


24.715 


16.915 


358.587 


7602 


7291 


10 


2.334 


8.375 


30.049 


26.382 


18.791 


426.442 


9677 


9391 



Table 2. The values for NLO and NNLO QCD contributions used in our fits and the 
N^LO Fade estimates. 



In the same Table, we also give the estimates for r{n) and C^^-* (n) , obtained by using the 
expanded [1/1] and [0/2] Fade approximants formulae of Eqs.(16)-(19). For completeness, 
in the last two columns of Table 1 we also present the estimates for N^LO contributions 
to the anomalous dimension function ^j!^g{As), obtained with the help of the expanded 
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[1/1] and [0/2] Pade approximants. One can see that the results of appUcations of [1/1] 

(3) 

and [0/2] Pade approximants for 7]v5(^) a-rs almost identical to each other. 

Using the numbers presented in Table 1, one can construct Pade motivated expressions 
for r(n) by substituting the estimates for Jj^gin) into Eq.(ll). It should be stressed, that 
the obtained estimates for r(n) qualitatively agree with the ones presented in Table 2 
within the "Pade world" only, namely only in the case of application, in Eq.(ll), of the [1/1] 
or [0/2] Pade estimate for the four-loop coefficient of the QCD /3-function P^. However, 
in the case of / = 4 the direct application of the [1/1] and [0/2] Pade approximants 
underestimates the calculated value of fS^ by a factor of over 2.5 (/33|[i/i] ~ 3217; /?3|[o/2] ~ 

3058). In view of this, the application of Eq.(ll) with the Pade estimated values of Jj^sin) 
and the explicit expression for /33-coefficient give estimates of r{n), drastically different 
from the ones presented in Table 2 (for example, for the case of application of [0/2] Pade 
estimates it gives r(2) 16.6,....,r(10) ^ 49.3). 

It is already known that the accuracy of estimates of the N'^LO coefficient of the 
QCD /3-function can be improved by some additional fits of the polynomial dependence 
of on the number of flavours / and by applying the asymptotic Pade approximant 
(APAP) formula [38|. Therefore, it might be interesting to consider the possibility of 
making Pade estimates of N^LO contributions to 7Ars(^s) (see Table 1) more theoretically 
motivated. Analogous steps were already done in Ref.[|^] in the analysis of the status of 
N^LO Pade estimates for the anomalous dimension function of quark mass. The agreement 
of the obtained estimates with the calculated four-loop QCD results of Ref.|40, 41 1 turned 

out to be reasonable. One can hope that the application of similar procedure for the 

(3) 

APAP estimates of 7]Y^-terms and the substitution of the results obtained in Eq.(ll) 
(together with the explicit expression for the /33-term) might improve the agreement with 
the estimates presented in Table 2. At this step we consider the estimates presented in 
Table 2 as suitable results for modelling the unknown effects of the N^LO corrctions, which 
depend on the N^LO expression for the coupling constant As. 

Within this approach, the uncertainties of the results of NNLO fits can be estimated by 
modelling q{n) and C^'^\n) by using the [0/1] Pade approximants, which give (7(n)|[o/i] = 
[p(n)]^ and C'*-^H'^)l[o/i] = [C^^\n)]'^. The estimated values of g(n)|[o/i] are correct in sign 
for n > 2, while for C^'^\n) the same feature takes place in the case of n > 4 moments. 
Moreover, for n > 4 the relative values of ratios 9(?T')|[o/i]/'7(^) are varying from 1.3 to 1.5, 
while similar ratios for NLO contributions to the coefficient function C(^^(n)|[o/i]/C'(^^(n) 
are changing from 1.6 at n = 4 to 1.2 at n = 10. This precision seems to be rather 
acceptable for the [0/1] Pade estimates, which in the case of each concrete fixed value of 
n are based on one input term of the corresponding perturbative series. 

It should also be stressed that the uncertainties in values of r(n) are not so important, 
since the results of our fits are more sensitive to the form of predictions of the Pade 
approximations for the N^LO contributions to the coefficient function (namely, C^^\n)- 
terms) . 

Prom the results presented in Table 2 one can conclude that the theoretical series 
for for large n {n > 4), relevant to the behavior of xF^^XjQ^) in the intermediate 
and large x-region, probably have sign constant structure with asymptotically increasing 
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positive coefficients. Therefore, the apphcations of the expanded [1/1] and [0/2] Pade 
approximants for estimating the terms C^^\n) with n > 4 (which in both cases have the 
same positive sign and the same order of magnitude) might be considered as the useful 
ingredient for the N^LO Pade-motivated fits. 

However, in the cases of coefficient functions of n = 2, 3 moments of 2:^3 our intuition 
suggest nothing about the sign and order of magnitude of the third term in perturbative 
series cj^^(^) = 1 - 1.78^, - 47.47^2 ^nd cj^^(^) = 1 + 1.67A, - 12.71^^. Indeed, in 
these two cases the manipulations with [1/1] and [0/2] Pade approximants give drastically 
different estimates for the terms C^^\n) which for n = 2,3 differ both in sign and size 
(see Table 2). It is possible that this feature is related to the fact that for n = 2,3 
the coefficients of C^^{As) have no immediate (+l)'"m! growth, but exhibit some zigzag 
structure manifesting itself in the cases of definite perturbative series of quantum field 



theory models (for discussions, see e.g. Ref.[42|). This might give additional theoretical 
uncertainties of modelling higher-order perturbative QCD predictions for xF3{x,Q'^) in 
the region of relatively small x. 

In view of the questionable asymptotic behavior of the NNLO series for coefficient 
functions of the NS moments with low n (n = 2,3), we also use the idea of Ref.[35] and 



consider non-expanded Pade approximants in the process of analysis of the DIS data. 

Let us recall that the corresponding non-expanded [1/1] Pade approximants can be 
defined as 

i -|- Ag 

^i.1(^^)l[i/i]=f^4^ (21) 
l + d\ 'As 

whereaS") = (^[p(n)]2-g(n)^ /p(n), 6^"^ = -q(n)/p(n) and = (^[C«(n)]2-C7(2)(n)^/C«(n) 

4") = _C(2)(n)/C«(n). 

The explicit expressions for the non-expanded [0/2] Pade approximants read: 

^D(n,^)|ro/2i = n — (22) 



CiVs(^) I [0/2] = , . ,(n) . , ,(n) (2^) 



^''^''^■^"'"'^'"i + <M. + 4"^A^ 

where 6^"^ = -p{n), 6^"^ = p{nf - q{n), = -C^^\n) and df = [C^'^\n)f - C(2)(n). 
Since we consider the applications of both [1/1] and [0/2] Pade approximants as attempts 
to model the behavior of the perturbative series for the NS Mellin moments beyond the 
NNLO level, we use, in Eqs.(20)-(23), the N^LO expression for the coupling constant 
As, defined through Eqs.(12)-(14). It is worth mentioning here that quite recently the 
expanded and non-expanded Pade approximants were successfully used to study the N^LO 
approximation of the ground state energy in quantum mechanics [B3] and the behavior of 



the /?- function for the quartic Higgs coupling in the Standard Electroweak Model [44] 
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The next step is the reconstruction of the structure function xF3{x,Q'^) with both 
target mass corrections and twist-4 terms taken into account. The reconstructed SF can 
be expressed as: 

^ max 

xFg^— (x,Q2) = w{a,f3){x)Y.Q'^J{x)Y,cf\a,f3)M3+2,.F,{Q^) (24) 

n=0 j=0 

Mxi 

where Q^'^ the Jacobi polynomials, c~^\a,l3) contain a and j3 dependent Euler F- 
functions where a, (3 are the Jacobi polynomials parameters, fixed by the minimization of 
the error in the reconstruction of the SF, and w{a,[3) = x"(l — x)^ is the correspond- 
ing weight function. To take into account the target mass corrections, the Nachtmann 
moments 

Mn,.F. - M^^gm = I" ^^3(x,Q^)i±^^, (25) 

can be used, where ^ = 2x/(l + V), V = + 'iM'^^^x'^ / Q"^ and Mnud is the mass of a 
nucleon. However, to simplify the analysis, it is convenient to expand equation (|2^) into a 
series in powers of M'^^^JQ'^ [^. Taking into account the order 0{M^^^JQ'^) corrections, 
we get 



2\ 



(„ + 2)(n + l)n M^^,, .^2^ , ^.^n«d ^ 



2(n + 4) Q4 n+my-^ ^ ' q, 

We have checked that the order 0{M^^^i/Q^) terms in Eq.(26) have a rather small 
effect in the process of concrete fits. Therefore, in what follows we will use only the first 
two terms in the r.h.s. of Eq.(26). 

The form of twist-4 contributions h{x) in Eq.(24) was first fixed as 



Q 

where 



7 / \ Nmax (^) 



3 

n=0 j=0 



Mi%{Q') = C(n)M„^,%(Q2)^ + O(i^) (28) 



with A2 taken as a free parameter and C{n) defined following the IRR model estimates 
of Ref.[|l9| as C{n) = -n - 4 + 2/(n + 1) + 4/(n + 2) + 45i(n) (5i(n) = ^2]=! Vj)- It 
should be stressed that the multiplicative QCD expression M^^p^{Q'^) in Eq.(28), gener- 
ally speaking different from the intrinsic coefficient function of the twist-4 contribution. 
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leads to theoretical uncertainties in the contributions of higher-order QCD corrections to 
the twist-4 part of xF3{x,Q'^). This could provide additional theoretical errors in the 
studies of the status of the IRR-model predictions for the twist-4 terms in NNLO. 

To analyze this question at a more definite theoretical level, it is instructive to model 
the function h{x) by additional free parameters of the fits, not related to the IRR-model 
predictions. 

We will estimate the uncertainties of the values of A^^, and thus as{Mz), by studying 
the factorization and renormalization scale dependence of the outcomes of the fits. We will 
also analyze the stability of extracted values of the IRR-model parameter A2 and the twist- 
4 function h{x) to the explicitly calculated N^LO QCD corrections and other unknown 
N^LO terms (modelled with the help of the Fade resummation technique) included into 
the fits of concrete data. Our aim will also be the study of the influence of the choice 
of the initial scale on the results of Ref.|ll| and especially on those, which describe the 
x-shape of h{x) for xF^ within the method adopted by us. 



3 (a). The analysis of the experimental data: the extraction of A^g^ vs 



value. 

The results for our NLO and NNLO fits, made for the case of number of active flavours 



/ = 4, are presented in Table 1 of Ref.|ll|, where the values of the parameters for the model 
of xF^ ^, 6, c, 7 7^ (related to the parton distribution parameters) are also given. The 
results of Ref.0 were obtained by using the fixed value of the initial point Qq = 5 GeV^ 
and the fixed weight function of the Jacobi polynomials reconstruction formula of Eq.(24), 
namely x^'"^ {1 — x)^ . Note, that this form is similar to the x-shape of the NS structure 
function itself. Indeed, the value of the parameters /3 = 3 is in agreement with the quark- 
counting rules of Ref . p6[ , while the value a = 0.7 is close to the value of the parameter, 
which describe the Regge theory behaviour of NS SF at small x. 

In this section, we will present more definite arguments in favour of the used form of 
the Jacobi polynomial weight function and will study the initial Q'q- scale dependence of 
the results for A^^, extracted in different orders of perturbation theory. 

We will also construct the N^LO Q^-evolution equations for the Mellin moments using 
the Fade approximants, written down both in the expanded and non-expanded forms. 
In the process of these "approximate" N^LO fits, the explicit N^LO expression for the 
QCD running coupling constant As, defined in Eqs.(12)-(15), will be used. Thus, from 
the fits we obtain the N^LO estimates of the parameter A^^ (and therefore as(M^)), 

and of the common factor ^2 of the IRR model. The comments on attempts to apply the 
scheme-invariant analysis for estimating the renormalization-scheme dependence of the 
results obtained will also be presented. 

It should be stressed that despite the general theoretical preference of applications of 
the diagonal Fade approximants (for a recent analysis, see, e.g., Ref.[^), the N^LO [1/1] 
Fade approximant description of the CCFR'97 experimental data turned out to be not 
acceptable in our case, since it produces a rather high value of IK^ /nep > 2 (where 
nep = 86 is the number of experimental points, taken into account in the case of the cut 

> 5 GeV"^). However, the application of [0/2] Fade approximants produced reasonable 
results. We think that the non-applicability of the [1/1] Fade method in the process of 
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fitting CCFR xF^ data with the help of the Jacobi polynomial approach can be related to 
a rather large value of the ratio [C(^'*(2)/C*^^^(2)]|[x/i] in the expression for the NS moment 
Mi^xF3 (see Table 2). 

A similar effect of the preference of the [0/2] Fade approximant analysis over the [1/1] 
one was found in Ref. p^ from the comparison of QCD theoretical predictions for the 
polarized Bjorken sum rule (which are closely related to the QCD predictions for the 
first moment of 2:^3, namely, for the Gross-Llewellyn Smith sum rule) with the available 
experimental data. 

Considering the problem of minimization of the dependence of the results of the fits 
on free parameters a, (5, we found several minima on the (a, /?)- plane at Qo = 5 GeV^: 



1. Minimum A: a//3 ss-0.6/0.55. 

At this minimum we got reasonable LO, NLO and NNLO values of A-^i^^ for N^ax = 
6. However, the appearance of this minimum strongly depends on the number of 
moments taken into account. For example, in the case of N^ax = 10 we were unable 
to find this minimum at LO and NLO , so we consider this minimum as a spurious 
one; 

2. Minimum B : a//3 ps-0.5/-0.9. 

This minimum appears in LO and NLO. However, this minimum does not appear 
in NNLO, so we consider it as non-applicable for our NNLO fits. 

3. Minimum C: q//3 ssO. 8/1.3. 

For Njnax = 6 the LO and NLO values of can be obtained. However, in NNLO 
this minimum does not manifest itself. Moreover, it disappears in LO and NLO for 
the case of N^ax = 10. Therefore, we consider it as the spurious one also. 

4. Minimum B: a/(3^ 0.6/-0.99. 

It should be stressed that in LO and NLO this minimum appears only for N^ax = 6 
and disappears for Nmax = 10. Moreover, the obtained value of (5 results in the 
unnatural singular 1/(1 — x) behaviour of the Jacobi polynomial weight function 
w{a,f3). In view of this we consider this minimum as the unphysical one. 

5. Minimum E: a//5 ~ 0.7/3.0. 

It is the minimum at which we worked earlier in Refs.[pT|, |l2|. At this minimum the 



results of LO and NLO fits are in agreement with the ones obtained in Refs.[ll, 12|. 
It should be noted that in LO and NLO this minimum is stable due to variation of 
Qq and to the inclusion of higher Mellin moments into the reconstruction formulae of 
Eq.(24) and Eq.(27) (we checked this statement, repeating the fits for Nmax = 10). 
In NNLO this minimum appears at higher than 10 GeV^. 

Since the values of parameters a, /? at Minimum E are identical to the initially consid- 
ered ones( a = 0.7, /3 = 3) and in view of the stability of the results of the LO and NLO 
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fits to the value of the initial scale Q'q, we consider this minimum as the physical one and 
will work in its vicinity , fixing a = 0.7 and /3 = 3 as earlier. 

To study the dependence from the choice of the initial scale in more detail and thus 
to check the reliability of the results, obtained in Refs.[0, 12|, we performed LO, NLO 



and NNLO fits for different values of the initial scale without taking account of twist-4 
contributions, but with target mass corrections included (see Ref.|0). The results for 
A^IL are presented in Table 3. 

MS ^ 





5 


8 


10 


20 


50 


100 


LO 


266±35 


266±35 


265±34 


264±35 


264±36 


263±36 


LO* 


382±38 


380±41 


380±40 


379±46 


378±43 


377±42 


NLO 


341±30 


340±40 


340±35 


339±36 


337±34 


337±37 


NLO* 


322±29 


321±33 


321±33 


320±34 


319±36 


318±36 


NNLO 


293±30 


312±33 


318±33 


326±35 


326±36 


325±36 


NNLO* 


284±28 


312±33 


318±33 


326±35 


326±36 


325±36 


N^LO [0/2] 


293±29 


323±32 


330±35 


335±37 


326± 36 


319±35 



Table 3. The Ql dependence of [MeV]. The LO* means that in the LO-fits NLO 
as is used; NLO* (NNLO*) indicates that in the NLO (NNLO) fits NNLO (N^LO) is 
used. The N^LO [0/2] marks the results of the expanded [0/2] Fade fits with defined 
in N^LO. 



One can see that the LO and NLO results are stable to the variation of Qg. The 
results of LO* fits are higher than the LO ones, and from this level other perturbative 
QCD effects tend to decrease the values of Iv^YTa to the level of the NNLO*-fits. 

M S 

The NNLO results are sensitive to the variation of the initial scale Qq. The values of 
AgL become stable for Ql > 10 GeV"^ only. The same effect manifests itself for the results 
of the expanded [0/2] Fade approximant fits, which incorporate the explicit information 
about the N^LO expression for the coupling constant Ag. We think that this effect might 
be related to a rather peculiar behaviour of the NNLO perturbative QCD expression of 
n = 2 moment. Indeed, taking into account the exact numerical values of the coefficients 
p(2), q{2), C^^\2) and C'-^^(2) from Table 2, we find that the perturbative behaviour of 
this moment is determined by the following perturbative series 

AD{2,Ag)d^l{As) = 1 - 0.132^(Q2) _ 46.155^^(^2) ^ _^ ^29) 

where the relatively large coefficient is dominated by the NNLO term of the coefficient 
function of n = 2 moment. Thus we think that it is more appropriate to start the NNLO 
QCD evolution from the initial scale Ql = 20 GeV^, where the numerical value of the 
contribution in Eq.(29) is smaller. Note that this choice of the initial scale is also 
empirically supported by the fact that it coincides with the middle of Q^ kinematic region 
of the CCFR data. 

In Table 4 we present the results of our new fits for and IRR-model parameter A'2 
obtained in LO, NLO, NNLO and N'^LO ( modelled by the expanded and non-expanded 
Fade approximants) in the cases of both 7 7^0 and 7 = 0. 
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7 - free 


7 = 0- fixed 


Q2 > 


^% (MeV) 


A^(HT) 


/points 


(MeV) 


^^(HT) (GeV^)) 


/points 


5 GeV'^ 














LO 


264±35 




113.1/86 


241±35 




121.7/86 




433±52 


-0.33±0.06 


82.8/86 


398±71 


-0.31±0.08 


121.7/86 


NLO 


339±36 




87.6/86 


313±36 




95.2/86 




369±45 


-0.12±0.06 


82.3/86 


341±36 


-0.11±0.05 


92.3/86 


NNLO 


326±35 




77.0/86 


314±35 




86.1/86 




326±35 


-0.01±0.05 


76.9/86 


315±34 


-0.02±0.05 


86.3/86 


N^LO 


332±28 




76.9/86 


314±28 




86.3/86 


(n.e.) 


333±27 


-0.04±0.05 


76.3/86 


315±27 


-0.04±0.05 


85.7/86 


N^LO 


335±37 




77.9/86 


328±37 




85.1/86 




340±37 


-0.04±0.05 


77.2/86 


335±37 


-0.05±0.05 


84.2/86 


10 GeV^ 














LO 


287±42 




77.3/63 


283±39 




78.1/63 




529±77 


-0.52±0.12 


57.8/63 


515±75 


-0.50±0.12 


59.5/63 


NLO 


349±40 




63.9/63 


344±44 




64.8/63 




436±55 


-0.24±0.10 


58.3/63 


427±55 


-0.23±0.10 


59.5/63 


NNLO 


338±30 




57.4/63 


337±40 




58.7/63 




354±45 


-0.06±0.09 


56.9/63 


352±42 


-0.03±0.09 


58.2/63 


N^LO 


348±41 




57.3/63 


347±41 




58.0/63 


(n.e.) 


373±48 


-0.09±0.09 


56.2/63 


373±48 


-0.09±0.09 


56.9/63 


N^LO 


344±40 




56.8/63 


345±30 




57.3/63 




3(j2±iG 


-0.07±0.(J9 


56.2/63 


363±-lG 


-0.07±0.09 


56.6/63 


15 GeV'^ 














LO 


319±48 




58.5/50 


320±47 




58.5/50 




530±89 


-0.56±0.18 


49.9/50 


525±45 


-0.56±0.20 


50.0/50 


NLO 


365±46 




52.3/50 


366±46 




52.3/50 




440±71 


-0.25±0.17 


50.3/50 


438±69 


-0.24±0.17 


50.3/50 


NNLO 


343±44 




50.4/50 


341±44 




50.8/50 




350±56 


-0.03±0.15 


50.3/50 


347±55 


-0.03±0.15 


50.8/50 


N=^LO 


348±45 




50.3/50 


348±45 




50.4/50 


(n.e.) 


358±57 


-0.04±0.14 


50.2/50 


357±56 


-0.04±0.15 


50.3/50 


N^LO 


342±43 




50.2/50 


343±41 




50.3/50 




347±58 


-0.02±0.15 


50.2/50 


347±54 


-0.02±0.15 


50.2/50 



Table 4. The results of extractions of the parameter A-^^ and the IRR coefficient A'2, 
(in GeV^) defined in Eq.(23), from LO, NLO, NNLO and N^^LO non-expanded (n.e.) and 
expanded Fade fits of CCFR'97 data. In the fits we take Qg=20 GeV^ 
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Looking carefully at Table 4 we arrive at the following conclusions: 



• The results of LO and NLO fits are identical to the ones obtained in Ref. . 

• Our fits demonstrate that the NNLO values of the parameter depend on the 
choice of the initial scale Q'q. In the case of Q'q = 20 GeV^ the NNLO perturba- 
tive QCD contributions are less important, than in the case Qq = 5 GeV^, earlier 
considered in Ref. . Indeed, for different Q^-cuts they change slightly the NLO 
values of A.^^, provided twist-4 corrections are switched off. 

• As was mentioned previously, this effect might be related to the peculiar behaviour 
of the NNLO perturbative expression for n = 2 moment of xF^ (see Eq.(29)) and, 
therefore, to the theoretical uncertainty of the NNLO behaviour of xF^ at small 
X. We checked this conclusion by comparing the results of the NLO and NNLO 
Ql = 5 GeV'^-fits of the CCFR'97 data cut at x > 0.04. The result of these test-fits 
demonstrate the tendency, identical to the one revealed after moving the initial scale 
to Qo — 20 GeV^ , namely, the minimization of the difference between the values of 
A^^, extracted in NLO and NNLO. 

• However, in the case when the IRR- model for the twist-4 corrections are included 
into the analysis, the effects of the NNLO corrections are still important and decrease 
both the value of A^^ and the IRR- model parameter A'2, making the first one almost 

identical to the NNLO value of A^^, obtained without twist-4 terms. The latter 

MS' 

one is compatible with zero within statistical error bars. A similar feature was also 
observed in the case of the fits, made in Ref.|pT| for the initial scale Qq = 5 GeV^. 
This property confirms the conclusion of Ref. that the results of the NNLO fits 
are less sensitive to the parameter of the IRR-model of the twist-4 term. A similar 
conclusion was also recently made in Ref. ||4^ while comparing the experimental 
data for the DIS R-ratio with the available NNLO perturbative QCD results of 
Ref.|49], although the earlier analysis of the experimental data for this quantity 



with a different kind of the HT model leaves still room for the power suppressed 
behaviour [Eol, 51]. 



The values of A^i^, obtained from the fits with the [0/2] Fade estimates (both 
in the expanded and non-expanded variants) turn out to be almost insensitive to 
the choice of the Q^-cut of the data, the value of 7 and thus incorporation of the 
(l+7x)-factor in the parton distribution model. The latter fact, in turn, can indicate 
that the change of the model xF^ix, Qg) = A{Ql)x^'^Qo\l - x)<^o\l + -/{Ql))x) to 
xF3ix,Ql) = A{Ql)x^^Qo){l - x)<Qo){l + j{Ql)x + e{Ql)^), used in the MRST 
and GRV fits, might affect the obtained results only slightly; 

Large errors in definite results for A^i^, presented in Table 4, reflect the correlations 
of these uncertainties with the errors of the IRR-model parameter A'2; 
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• The property Xlo > Xnlo > Xnnlo ~ x'n^lo I'^flects the importance of the effects 
of higher order perturbative QCD corrections in the process of fits of the concrete 
experimental data; 

• For ah Q^-cuts the expanded N^LO results for A^l^ are almost identical to the ones 
obtained with the expansion of the Fade approximants in Taylor series. Moreover, 
the x^-criterion does not discriminate between these two variants of the Fade mo- 
tivated fits (see especially the results, obtained for the cuts > 10 GeV^ and 

> 15 GeV'^). In our future studies we will consider the results of applications of 
the expanded Fade approximants. 

The results of the the NNLO fits, made with the help of the Jacobi polynomial ex- 
pansions method, are compared to the available CCFR'97 data in Fig.l. Drawing the 
theoretical curves we used Eq.(24) with zero twist-4 contributions. The value of the QCD 
scale parameter, which governs the theoretical behaviour of the moments of xF^ SF, turned 
out to be A^i^ = 326 ± 35 MeV (the of the corresponding fits is 77.0/86). The values 
of the corresponding parameters of the xF^ model at Qq = 20 GeV^ are : A = 4.70 ±0.34, 
b = 0.65 ± 0.03, c = 3.88 ± 0.08, 7 = 0.80 ± 0.28. One can see that the NNLO results of 
the fits without twist-4 corrections are in good agreement with the CCFR'97 experimental 
data for xF^, . 

To determine now the values of as{Mz), we transformed through the threshold 
of the production of the fifth flavour, M5. This is done using the LO, NLO, NNLO and 
N^LO variants of the MS'-scheme matching conditions, derived in Ref.|^^ following the 
lines of Ref . j23U . The related values of Apr^ can be obtained with the help of the following 
equation: 
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Fig.l The comparison of the CCFR'97 data with the results of our NNLO Jacobi 

polynomial fits. 
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where C2 = —7/24 was calculated in Ref.[^ (see also Erratum to Ref . p^] ) and the 
analytic expression for C3, namely C3 = -(80507/27648)C(3) - (2/3)C(2)((l/3)/n2 + 
1) - 58933/124416 + (//9)[C(2) + 2479/3456] was recently found in Ref.@. Here /?/ 
if^i~^^) coefficients of the /^-function with / (/ + 1) numbers of active flavours, Lh = 
/n(M|^-^/A^^^) and Mj+i is the threshold of the production of a quark of the (/ + 1) 
flavour. In our analysis we will take / = 4 and nib ~ 4.8 GeV and vary the threshold of 
the production of the fifth fiavour from M| = to M| = (6mf,)^ in accordance with the 
proposal of Ref. 1 26]. The difference between different prescriptions of fixing the matching 
point is included into the estimate of the theoretical uncertainties of the final results for 
as{Mz). 

In the case of the non-zero values of the twist-4 function h{x) 7^ 0, the results of the 
fits are presented in Table 5 in the next Section. 

It should be stressed that we consider the outcomes of our N^LO approximate fits as 
theoretical uncertainties of the NNLO results in the same manner as the results of the 
NNLO analysis are considered as the measure of theoretical uncertainties of the NLO 
results. In particular, we introduce the characteristic deviations l\^^^o_ |(^^W_)^^iO _ 

^a{4) \NNL0\ aNLO— I('a('') \NNL0 _ /a(4) \NL0\ 

When the twist-4 terms are included into the fits, the difference /\^^^0— | (^\^^_^n^lo _ 
{K^)NNLO\ ig smaller than the NLO correction term A^^^^ \{}M\NNL0 _^^W_^nlo\^ 
A similar tendency ^ ^nlo ^^^kes place in the case of the fits without twist-4 

corrections. These observed properties demonstrate the reduction of theoretical errors due 
to cutting the analyzed perturbation series in different orders. 

It is known that the inclusion of higher-order perturbative QCD corrections in the 
comparison with experimental data decreases the scale-scheme theoretical errors of the 



results for and thus as{Mz) (see e.g. Refs.|5^, |l^, 34]). Among the ways of prob- 



ing the scale-scheme uncertainties are the scheme-invariant methods, namely the principle 
of minimal sensitivity, the effective charge approach ( which is known to be identical to 
the scheme-invariant perturbation theory) and the BLM approach (for a review of these 
methods see e.g. Ref.[^). The scheme- invariant methods were already used to estimate 
the unknown higher order corrections in SFs (see Ref. 1 51], where a strong decrease of the 
value of the QCD scale parameter was found in the NLO scheme-invariant fit of the exper- 
imental data for the NS part of the F2 ). It was also used to try to estimate the unknown 
at present N^LO corrections to definite physical quantities |36], and DIS sum rules among 
others. Note that the predictions of Ref.p^ turned out to be in agreement with the re- 
sults of applications of the Fade resummation technique (see Ref.||3^). Therefore, we can 
conclude that the application of the methods of the Fade approximants should lead to 
the reduction of the scale-scheme dependence uncertainties of the values of as{Mz) in the 
analysis of the CCFR data. 

To consider the applicability of the Fade resummation technique for fixing scale-scheme 
dependence ambiguities, we performed the scheme-invariant fits following the ideas of 
Ref. []5^ . We found that in NNLO the application of the effective charge approach gave a 
rather high value of (x^ ~ 111/86). This, in a turn, can be related to the appearance of 
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large and positive values of the NNLO terms /?2(?T-)e// of the effective-charges /3-functions, 
which are the important ingredients of the scheme-invariant approach of Ref.||5j]. Similar 
problems have also been observed in the case of the scheme-invariant applications to 
the study of the NNLO perturbative QCD predictions to other renormalization-group 



invariant quantities (see Refs.[p5|, 41, |5^ for discussions). In this work we avoid the 
detailed investigation of this problem. However, in Sec. 4 we will discuss the results, which 
were obtained after variation of the factorization and renormalization scales without using 
the freedom in the choice of the scheme-dependent NLO and NNLO corrections to the 
anomalous dimension functions and NNLO corrections to the QCD /3-function. 

We now present the values of as{Mz), extracted from the fits of the CCFR'97 exper- 
imental data for xF-^, obtained with the twist-4 contribution, modelled through the IRR 



model of Ref.119 



NLO HT of i?e/.[19] as[Mz) = 0.120 ± 0.003(staf) (34) 

±0.005(syst) ± Q.QM{theory) 

NNLO HT of -Re/.[19] as{Mz) = 0.118 ± 0m2{stat) (35) 

±0.005(syst) ± 0.003(t/ieory) . 

Anticipating the considerations of Sec. 3(b) we also present the results of NLO and 
NNLO extractions of as{Mz) with the twist-4 contribution, modelled by additional free 
parameters of the fit: 

NLO HT free a,{Mz) = 0.123l[}:[]?g (stat) (36) 
±0.005(syst) ± 0.004(t/ieory) 

NNLO HT free as{Mz) = 0.12lt°:°°^(stat) (37) 
±0.005(syst) ± 0.003(t/ieory) . 

Systematic uncertainties are taken from the CCFR experimental analysis, presented in 
the first work of Ref.|0, and the theoretical uncertainties in the results of Eqs.(34),(36) 
[Eqs.(35),(37)] are estimated by differences between central values of the outcomes of the 
NNLO and NLO [N^LO and NNLO] fits, presented in Tables 4,5, plus in the apphcation 



of the MS-scheme matching condition (which following the considerations of Ref.|26] was 
estimated as ^as{Mz) = ±0.002). In the process of fixing the theoretical errors with 
the consideration of the N^LO corrections we take into account the differences between 
applications of the expanded and non-expanded Fade approximants. 

It can be seen that due to a large overall number of the fitted parameters the results 
of Eqs.(36),(37) for as{Mz) obtain large statistical uncertainties. As can be seen from the 
results of Eq.(34),(35) for the QCD coupling constant, it is possible to decrease their values 
by fixing the concrete form of the twist-4 parameter h{x). However, if one is interested 
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in the extraction of the form of the twist-4 parameter h{x), one should take for granted 
these intrinsic theoretical uncertainties of the value of as{Mz)- 

3 (b). The analysis of the experimental data: the extraction of the shape 
of the twist-4 terms. 

Apart from the perturbative QCD contributions, the expressions for DIS structure 
functions should contain power-suppressed high-twist terms, which reflect possible non- 
perturbative QCD effects. The studies of these terms have a rather long history. At 
the beginning of these studies it was realized that the twist-4 contributions to structure 
functions should have the pole-like behavior ~ 1/((1 — x)Q^) |5^, 55]. This behavior 



was used in the phenomenological investigations of the earlier less precise DIS uN data 
|60| , SI I , which together with other different procedures of analyzing neutrino DIS data 
[|62|, B^] was considered as the source of the information about scaling violation parameters. 



The development of the renormalon technique (See Refs.[21, 19, |20[] and Ref.[22| for a 



detailed review) pushed ahead a more detailed phenomenological analysis of the possibility 
of detecting higher- twist components in the most precise DIS data available at present, 
obtained by BCDMS, SLAC, CCFR and other collaborations. It turned out that despite 
the qualitative status of the renormalon approach, a satisfactory description of the results 
of the QCD NLO F2 analysis [^] in terms of the IRR technique was achieved [H, p. The 
next step was to clarify the status of the predictions of Ref . [p!S|] for the form and sign of 
the twist-4 contributions to xF^ . The study of this problem was done in Ref. [|ll| (see also 
Ref. [p!o[ ) . In this section, we discuss the results of a more refined analysis of the behavior 
of the twist-4 contributions to xF^ in LO, NLO, NNLO and beyond. 

In Table 4 we present the dependence of the extracted value of the parameter A2 on 
different orders of perturbative QCD predictions, Q^-cuts of the CCFR experimental data 
and the coefficient 7 of the parton distribution model for xF^, fixing the factorization 
scale Qq = 20 GeV^. Note that the parameter A2 was introduced in the IRR model of 
Eq.(28), taken from Ref.||], and fixed there as A'2 « -0.2 GeV^, which is necessary for 
the description of the fitted twist-4 results of Ref. 1 64] for F2 in the IRR language. We 
found that the value of this parameter, extracted in LO and NLO, is negative, differs from 
zero for about one standard deviation and qualitatively agrees with the IRR-motivated 



guess of Ref. 119 1 . Moreover, the results of our LO and NLO fits are also in agreement 
with the value of the parameter h = — 0.38 ± 0.06 GeV'^ of a different model of the twist-4 
contribution to xF^, namely xF^(x,Q'^)h/{{l — x)Q^), extracted previously in Ref. ]l6l|] 
from the old vN DIS data. 

It is interesting to notice that the results of Table 4 reveal that for larger Q^-cuts 
10 - 15 GeV^ the values of A'^ in the LO and NLO fits are less sensitive to the included 
number of experimental points than in the case of the low Q2 cut (5 GeV^). This feature 
can be related to the logarithmic increase of the QCD coupling constant Ag at lower 
Q^. However, since we are interested in the extraction of the power-suppressed twist-4 
contribution, we concentrate on the discussion of more informative, from our point of view, 
fits with low Q^-cut 5 GeV^, which contain more experimental points and thus are more 
statistically motivated. 

We now turn to a pure phenomenological extraction of the twist-4 contribution h{x) 
to xF^ (see Eq.(24)), which is motivated by the work of Ref. 1 64] for F2. In the framework 
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of this approach the x-shape of h{x) is parametrized by additional parameters hi = h{xi), 
where Xi are points of the experimental data bining. The results of the multiloop extrac- 
tions of these parameters are presented in Table 5 and are illustrated in Fig.2. 
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U.UoO 


-u.uio It u.zuo 


U.IOZ It u.zou 


0119-1-0 911 
U.llZ It W.Zll 


0110-1-01 Q'^ 

u.liu It u.iyo 


u.uou 


u.uoo It u.ioy 


U.ZOD It U.Z4tU 


n 91 /I -1- n 9nn 
U.Z14 It u.zuu 


1 71 -1- 1 RS 
U. 1 ( 1 It U. IDo 


U.U 1 u 


n 141 4- n 1 
-u.i^i in u.ioy 


U.l^^ It U.ZOO 


n 1 n« 4_ n 909 
U.±UO It u.zuz 


n 01 7 4- 1 c;i 
U.Ul 1 It U. 101 


0.090 


-0.177 it 0.127 


0.144 lb 0.254 


0.142 lb 0.202 


0.026 lb 0.142 


0.110 


-0.343 lb 0.127 


-0.009 lb 0.237 


0.045 lb 0.205 


-0.080 lb 0.146 


0.140 


-0.408 lb 0.116 


-0.085 lb 0.174 


0.060 lb 0.187 


-0.049 lb 0.140 


0.180 


-0.351 lb 0.173 


-0.072 lb 0.128 


0.154 lb 0.161 


0.093 lb 0.145 


0.225 


-0.547 lb 0.244 


-0.351 ± 0.180 


-0.098 ± 0.128 


-0.103 lb 0.127 


0.275 


-0.548 lb 0.334 


-0.472 lb 0.321 


-0.228 lb 0.185 


-0.193 lb 0.169 


0.350 


-0.295 lb 0.414 


-0.390 lb 0.470 


-0.154 lb 0.264 


-0.104 lb 0.234 


0.450 


-0.098 lb 0.410 


-0.307 lb 0.504 


-0.146 lb 0.330 


-0.125 lb 0.303 


0.550 


0.095 ± 0.324 


-0.134 lb 0.415 


-0.121 lb 0.323 


-0.138 ± 0.317 


0.650 


0.380 lb 0.211 


0.215 lb 0.262 


0.140 lb 0.245 


0.111 lb 0.251 



Table 5. The results of extractions of the HT contribution h{x) to xF^ and the parameters 
A,b,c,j with the corresponding statistical errors. The QCD fits of CCFR'97 data were 
performed taking into account TMC in LO, NLO {N^ax = 10), NNLO and N^LO {Nmax = 
6). In the latter case the expanded [0/2] Fade approximants were used. The fits are done 
for Ql = 20 GeV^. 
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Fig. 2. The x-shape of h{x) extracted from the fits of CCFR'97 data in the case of fixing 

the intial scale at Qq = 20 GeV'^ 

Looking carefully at Table 5 and Fig. 2 we observe the following features: 
1. The x-shape of the twist-4 parameter is not inconsistent with the expected rise of 



h{x) for X — > 1 58] in all orders of perturbation theory; 



2. the values of parameters h{xi) at the upper and lower points of the kinematic region 
(xi6=0.650 and xi =0.0125) are stable to the inclusion of higher order perturbative 
QCD corrections and application of the Fade resummation technique. At large values 
of X this feature is in agreement with the previous statement; 

3. the function h{x) seems to cross zero twice: at small x of order 0.03 and larger 
X about 0.4. It should be noted that the sign-alternating behavior of the twist-4 
contributions to DIS structure functions was qualitatively predicted in Ref . [|60|] ; 

4. in LO and NLO our results are in qualitative agreement with the IRR prediction of 
Ref.p9[| (for discussions, see Ref. 1 22]); 



5. in NNLO this agreement is not so obvious, though a certain indication to the mani- 
festation of the twist-4 term survives even in NNLO; 
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6. however, in NNLO we observe the minimization of the amphtude of the h{x)- 
variation. Thus we conclude that the inclusion of the NNLO corrections into the 
game might "shadow" the effects of the power suppressed terms in NNLO. This 
property was previously observed in LO as a result of the analysis of less precise 
DIS neutrino data in Ref.|59]. In the modern experimental situation, namely, in the 



analysis of more precise DIS neutrino data of the CCFR collaboration, we observe 
this feature in NNLO; 

7. we checked the reliability of this NNLO foundation by going beyond this perturbative 
approximation using the method of Fade approximants. The result of this analysis 
reveals a relative stability of the NNLO results for h{x); 

8. the property of minimization of the x-shape of h{x) in NNLO and N^LO is identical 
to the effect of decreasing the IRR model parameter A2 in NNLO and N^LO (see 
Table 4); 

9. these observed properties clarify why the results of the NNLO and N^LO fits for 

, presented in Tables 4,5 practically do not depend on the inclusion of the 
twist-4 contribution using the IRR model. Indeed, at this level the twist-4 terms 
manifest themselves less obviously. 

From our point of view, statements (8)-(9) reflect self-consistency of the results of our 
different fits with twist-4 terms included in different ways. 

4. The factorization and renormalization scales uncertainties. 

As it is known from the work of Ref. [|6^ ] it is rather instructive to study the sensitivity 
of the results to the variation of renormalization and factorization scales. We will study the 
question of factorization-renormalization scale dependence within the class of MS'-scheme 
only. This means that we will change only the scales without varying the scheme-dependent 
coefficients of anomalous dimensions beta-function. 

The arbitrary factorization scale is entering in the following equation: 



where fip is the factorization scale and 



l + kiAs{Qy^il)+k2Al{Q^/fil) 



(38) 



k2 = kf + -^ki 

Let us choose the factorization scale as fjij^ = m|j^^-F- 
Then we have: 



fei = -/3o/n(l^) (39) 

Pi 



ki = poln{kF) (40) 
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In this case after application of the renormahzation group equation and substitution 
of Eq.(38) into Eqs.(7,8) of Sec.2 we get 



exp 



where a = 7^^*^/2/30 and 



P{x) 



dx 



{As{Q''kF)T X AD{n,As{Q^kF)) 



(41) 



AD{n,As{Q^kF)) 



1 + 



p{n) + aki 



AsiQ^kp) 



(42) 



+ 



where 



/ N / N, / -.N A, o(a + 1) , 2 
q{n) + p{n)ki{a + 1) + ^kia + \ ' kj 

Po 



ki = ln{kp)Po 



AliQ^kp) 



(43) 



Now let us study the factorization and renormahzation scale dependence in the case 
when kji = kp = k (see, e.g., Ref . []6^] ) . In this case we should modify the coefficient 
function in Eq.(6) only as 



= 1 + c'^^\n)As{Q^k) + C(2)(n) + C^^\n) (3oln{k) 



(44) 



The commonly accepted practice is vary k in the interval 1/4 < < 4 (see, e.g., 
Ref.H) 

We repeated our fits, described in Sec. 3a, both without and with IRR model of the 
twist-4 terms in the cases of /c = 1/4 and k = 4. The obtained results are presented in 
Table 6: 



Order 


k 


Afc 


A'2 (HT) 


X'^/points 


NLO 


4 


116 




99.1/86 




4 


213 


- 0.22±0.06 


84.2/86 




1/4 


-61 




80.4/86 




1/4 


-99 


+ 0.02±0.05 


80.2/86 


NNLO 


4 


35 




83.5/86 




4 


66 


- 0.11±0.06 


83.5/86 




1/4 


-51 




87.3/86 




1/4 


-45 


+ 0.09±0.05 


84.5/86 



Table 6. The results of NLO and NNLO fits of CCFR'97 xF^ data for > 5 GeV"^ with 
different values of factorization/renormalization scales. Ak{MeV) = ^]^(^) ~ ^'Ms^^ ~ 
1). The value of the IRR model coefficient is given in GeV^. The initial scale is fixed as 
Ql = 20 Gey 2. 
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The factorization/renormalization scale ambiguities of the NLO (NNLO) values of 
are giving the estimate of theoretical errors of as{Mz), which can be compared 
with the ones, given in Sec. 3 (see Eqs.(34),(35)). Transforming the results of Table 6 into 
the related errors of as{Mz) with the help of NLO and NNLO variants of the matching 
condition of Eq.(30) we get: 

Aas{Mz)NLO = i^o? (45) 

(]\/r \ +0.003 
l\as[Mz)NNLO = -0.002 

in the case when HT are included through the IRR model and 



Aa,{Mz)NLO = t'oml (46) 

A^. (i\/r \ — +0.003 

l\as{Mz)NNLO — -0.002 

in the case when HT-terms are neglected. 

One can see, that the inclusion of the HT terms even through the definite model leads 
to the increase of theoretical uncertainties of the NLO value of the as{Mz)- However, 
for the value of /c = 1/4 the fitted value of A'2{HT) is lying closer to zero. This feature 
demonstrates once more that the inclusion of the NNLO corrections, modeled using the 
renormalization/factorization scale dependence have the tendency to decrease the fitted 
value of IRR model parameter. Another interesting observation is that at the NNLO the 
scale-dependence is drastically smaller. 

It should be stressed, that the results obtained without HT terms are in agreement 
with the results of the NLO and NNLO NS analysis of theoretical predictions for the NS 
part of F2 SF, which were obtained recently in Ref.p5|] with the help of DGLAP approach. 
Moreover, considering the case when kp 7^ kn and varying both parameters within the 
interval 1/4 < /ci? < 4 and 1/4 < /c/j < 4, we found that in NNLO the results of Table 
6 depend only slightly on the choice of the factorization scale and are mostly related to 
the different values of renormalization scale. This feature confirms our findings that the 
NNLO corrections to the coefficient functions are more important than the corrections of 
the same order to the anomalous dimensions functions. 

The next wellcomed feature is that the renormalization/factorization- scale estimate 
of theoretical uncertainty of the final NNLO results are in agreement with the theoretical 
uncertainty of the NNLO result of Eq.(35), estimated by means of application of the Fade 
approximation technique at the N^LO. This fact gives additional theoretical support in 
favour of the reliability of the NNLO value of as(M^), given in Eq.(35), which is one of 
the main results of our considerations. 

Let us discuss the renormalization/factorization-scale dependence of the results of the 
fits, performed in the case when HT terms were included as the free parameters and 7 = 
(for simplicity) . In this case we found the following scale-dependence of the central values 
of as{Mz), given in Eqs.(36),(37): 

Aas{Mz)NLO = tl'^ (47) 
Aas{Mz)NNLO = ±0.003 
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Fig. 3. The x-shape of h{x) extracted from the fits of CCFR'97 data in the case of fixing 
the initial scale at Qq = 20 GeV^ and changing factorization/renormalization scales. 

These errors demostrate the same feautures, as in the case of the inclusion of the 
HT-terms through the IRR-model, namely 

1) the reduction of theoretical errors in NNLO; 

2) the correspondence of the theoretical error in NNLO to the one, previously estimated 
with the help of simulation of the effects of the N^LO corrections by means of the Fade 
resummation technique; 

3) the agreement of the results of Eq.(47) with the independent estimates of factoriza- 
tion/renormalization scale uncertainties, recently obtained in the NS case in Ref . [|65|] . 

In Fig. 3 we present the investigation of the scale-dependence of the x-shape of the 
model-independent twist-4 coefficient h{x), which was obtained in Sec. 3(b) in LO, NLO, 
NNLO and Fade-motivated N^LO in the cases oi kp = = 1 (see Fig. 2). More definitely. 
Fig. 3 demonstrate, how the NLO and NNLO plots of Fig. 2 are changing in two cases, 
namely when kp = k^ = 1/4 and kn = kp = 4. 

These plots demonstrate us the following typical feautures of scale-dependence of the 
x-shape of the twist-4 contributions: 

• In the case oi kn = kp = A the result of the NLO analysis is closer to the LO 
behaviour of h{x), obtained in Fig. 2 in the case of kpt = kp = 1, while in the case 
oi kn = kp = 1/4 the NLO x-shape of the twist-4 correction is almost identical 
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to the one, obtained in Sec. 3(b) in NNLO and indicate the partial reduction of the 
twist-4 model-independent contributions. This means, that considering the renor- 
malizatin/factorization scale ambiguities we can simulate in part the effects of the 
NNLO contributions. This feature confirmes the findings of Ref.[^6|, obtained at the 
NLO within DGLAP approach. 

• In the case of kfi = kp = 1/4 both NLO and NNLO x-shapes of h{x) are in agreement 
with the results of applications of the N^LO Pade approximation technique, shown 
in Fig. 2. 

• In the case of kji = kp = 4 the NNLO x-shape of h{x) is also comparable with min- 
imization of the contribution of the twist-4 corrections. However, for these values of 
renormalziation/factorization scales the detailed structure of h{x) is less vivid and 
has large errors. This property is explained by the effective enhancement of the con- 
tributions of the NNLO corrections to the coefficient functions of the corresponding 
Mellin moments (see Table 2 and Eq.(44)) and large correlations of HT effets with 
the values of the parameter ^j^. 

5. The attempt of the inclusion of nuclear corrections. 

The effects of nuclear corrections are the remaining important source of uncertainties 
of the analysis of the DIS data. This is especially important for the experiments on heavy 
targets and in the case of CCFR data-on iron ^^Fe. 

The attempts to study these effects were done in Ref. []67| in the framework of the 
Deuteron-motivated model. A satisfactory QCD description of the CCFR data for xF^ 
was achieved due to the fact that in this case the nuclear effects do not exceed 5 % effect. 
However, a more realistic description of the nuclear effect for xF^ in the case of ^^Fe- 



target [68| revealed the appearance of new 1/Q and 1/M corrections to NS moments. 



which have the following form 

+ <Ap2>5p2M„^^(Q2) 

+^^|^n(n + l)M„^/,(Q^) (48) 

where M is the mass of a nucleon, and for ^^Fe the parameters of the nuclear model 
adopted in Ref.[||] are < e >^ -56 MeV, < > /(2M) « 35 MeV, < Ap^ >Fe^ 
—0.17 GeV and the derivative dp2Mn{Q'^) takes into account that the target momentum 
p can be generally off-mass-shell. This effect results in the following contribution |^] 

(49) 



which is independent of the nuclear content of the target. The numerical values of 5^2 
were also presented in Ref . [|6^] . 
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Note that the effects of nuclear corrections in DIS were also recently studied in Ref. | 
in the case of xF^ and in Refs. [f7C|] , [^] in the case of F2 (for earlier related works, see e.g., 
Ref. [|72|| ) . However, in our studies we concentrate ourselves on the consideration of the 
results of Ref. II . 

We included the corrections of Eqs.(48)-(49) into our fits and observed an unacceptable 
increase of value. We think that this can be related to a possible asymptotic character 
of the 1/M-expansion in Eq.(48), since the third term in the brackets of the r.h.s. of 
Eq.(48) becomes comparable with the first term (which is equal to unity) for n ~ 8 used 
in our fits. Note that the moments with large n are important in the reconstruction of the 
behavior of xF^ as x — > 1. This observed feature necessitates the derivation of the explicit 
expression for M^(Q'^) which is not expanded in powers of 1/M-terms. It should be added 
that the problem of possible asymptotic nature of the power suppressed expansions was 
mentioned in the case of Ellis- Jaffe and Bjorken DIS sum rules in Ref. (7^. 

Another possibility to explain the non-convergence of our fits with the nuclear correc- 
tions of Eq.(48) might be related to the fact that the parton distribution model for the 
nuclear SF xF^ can be different from the canonical model used by us[^]. In any case, 
we think that the study of the problem of possible influence of heavy nuclear effects on 
the results of fits of xF^ data is still on the agenda. 

Conclusion 

In this work we presented the results of extractions of as{Mz) and twist-4 terms from 
the QCD analysis of the CCFR data taking into account definite QCD corrections at the 
NNLO and beyond. Within experimental and theoretical errors our results for as{Mz) 
are in agreement with other extractions of this fundamental parameter, including its world 
average value as{Mz) = 0.118 ± 0.005. 

Our estimate of the NNLO theoretical uncertainties is based on application of the [0/2] 
Fade approach at the N^LO level and the analysis of factorization and renormalization 
scale ambiguities. The uncertainties of our NNLO analysis can be decreased after explicit 
NNLO calculations of the NS Altarelli-Farisi kernel. It should be added, however, that the 
NLO results, obtained by us both for the x-shape of twist-4 corrections and for the as{Mz) 
value are in good agreement with the results of the NLO DGLAF analysis of the CCFR'97 
xFs and F2 SFs data |7|, which gives as{Mz) = 0.1222 ± 0.0048 (exp) ± 0.0040 (theor) 

As to the twist-4 term, we found that despite the qualitative agreement of our NLO 
results with the IRR model prediction, at the NNLO level its x-shape tends to decrease and 
is stable to the application of the [0/2] Fade motivated N'^LO analysis, which is supported 
by the NNLO fits with fixing factorization and renormalization scales as = kp = 1/4. 

This feature can be related to the fact that the analysis of the CCFR data cannot 
distinguish the twist-4 l/Q^ terms from the NNLO perturbative QCD approximations 
of the Mellin moments. This possible explanation is similar to the conclusions of the 
LO analysis of the old less precise neutrino DIS data, made by the authors of Ref.|59|. 
It is worth to remind that they were unable to distinguish between LO logarithmic and 
l/Q^-behavior of the QCD contributions to Mellin moments of xF^. The experimental 
precision achieved in our days might move this effect to NNLO. Another possibility is 
that the inclusion of NNLO perturbative QCD contributions makes the extraction of the 
l/Q^-corrections within IRR model approach and by the model-independent way more 
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problematic (for discussions of the perturbatively based alternative of the IRR language 
within quantum mechanics model see Ref.[76|). 

Another related explanation of the decrease of the twist-4 terms in NNLO come from 
the partial summation of the definite terms of the asymptotic perturbative QCD series 
and that the increase of the order of perturbative QCD analysis effectively suppresses 
the remaining sum of the perturbative QCD contribution. One can hope that future 
experiments of the NuTeV collaboration will allow one to get new experimental data at 
the precision level, necessary for extracting a more detailed information about higher twist 
contributions. 

Note added 

After the technical part of this work was done, we learned about the work of Ref.|77| 
where the NNLO analysis of F2 SLAC, HERA and BCDMS data was performed both in 
the singlet and nonsinglet cases with the help of the method of Bernstein polynomials 
0. The main result of this work is the NNLO value as{Mz) = 0.1172 ± 0.0024, which 
is in agreement with our findings. In another recent work, the first steps towards the 
inclusion of NNLO corrections to the NS part of F2 and modelling the NNLO corrections 
to the kernel in the x-space were made [ 65 1 . Our analysis of factorization-renormalization 
scale uncertainties confirms their findings. It should be also noted, that quite recently the 



NNLO results of Ref. |16|| were confirmed [79| with the help of other methods. We hope 
that our possible future studies will allow us to generalize the Jacobi polynomial NNLO 
analysis presented in this work to the case of F2 SF. 
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